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Abstract. The spectra of very high energy 7-radiation from distant extragalactic objects suffer significant defor- 
mations during the passage of primary 7-rays through the intergalactic medium. The recently reported fluxes 
of diffuse infrared background radiation indicate that we detect, most probably, heavily absorbed TeV radiation 
from Mkn 421 and Mkn 501. This implies that the absorption-corrected spectrum of Mkn 501 may contain a 
sharp pile-up which contradicts to the predictions of the conventional models of TeV blazars, and thus may leads 
to the so-called "IR background- TeV gamma-ray crisis" . To overcome this difficulty, in this paper we propose two 
independent hypotheses assuming that (i) the TeV radiation from Mkn 501 has a secondary origin, i.e. it is formed 
during the development of electron-photon cascades in the intergalactic medium initiated by primary 7-rays; (ii) 
the pile-up in the source spectrum is a result of comptonization (in deep Klein-Nishina regime) of ambient op- 
tical radiation by an ultrarelativistic conical cold outflow (jet) with bulk motion Lorentz factor To > 3 x 10 7 . 
Within the uncertainties caused by the limited energy resolution of spectral measurements, the observed TeV ra- 
diation of Mkn 501 formally can be explained by the intergalactic cascade 7-rays, assuming however an extremely 
low intergalactic magnetic field in the direction to the source at the level of < 10 -18 G. We also demonstrate 
that the "bulk motion comptonization" scenario can quite naturally reproduce the unusual spectral features in 
the absorption-corrected TeV spectrum of Mkn 501, and briefly discuss the astrophysical implications of this 
hypothesis. 

Key words, galaxies: BL Lacertae objects: individual: Mkn 501 - cosmology: diffuse radiation - gamma rays: 
observations - gamma rays: theory 



1. Introduction 

The Diffuse Extragalactic Background Radiation 
(DEBRA) at infrared to ultraviolet wavelengths carries 
crucial cosmological information about the galaxy forma- 
tion epochs. It is believed that this radiation basically 
consists of two emission components produced by stars 
and partly absorbed and re-emitted by dust during the 
entire history of evolution of galaxies. Consequently, two 
distinct bumps in the spectrum of red-shifted radiation 
at near infrared (NIR) A ~ 1 — 2 /im and far infrared 
(FIR) A ~ 100 — 200 /im wavelengths, and a mid infrared 
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Direct measurements of the Cosmic Infrared 
Background (CIB) radiation contain large uncertainties 
because of heavy contamination caused by foregrounds 
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of different origin (for review see Hauser & Dwek 2001). 
Gamma-ray astronomy offers a complementary approach 
to derive information about CIB. Although this method 
requires certain model assumptions about the primary 
(un-absorbed) spectrum of 7-rays, it has an adequate 
potential for robust conclusions concerning the absolute 
flux and the spectrum of CIB. Moreover, the study of 
angular and spectral properties of 7-radiation from pre- 
dicted giant electron-positron halos surrounding powerful 
nonthermal extragalactic object s like AGN and radio- 
galaxies with known redshifts (Aharonian et al. 1994) 
can provide a unique tool to "measure" unambiguously 
the broad-band spectrum and the absolute flux of CIB 
at different cosmological epochs, and thus to probe the 
evolution of galaxies in past. 

The spectra of high energy (E > 10 GeV) 7-radiation 
observed from distant extragalactic objects suffer signif- 
icant deformation during their passage through the in- 
tergalactic m£dbrm_djjrj tojrrto 

with CIB (|Nikishov 1962] ; |Gould fc Schreder 1967| ; jSteckei 
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et al. 1992 ). The absorption features in the 7-ray spectra 
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depends on the flux of CIB, thus the study of such features 
from extragalactic objects with firmly established dis- 
tances could yield important constraints on CIB. Strictly 
speaking, this approach requires good understanding of 
the source spectra of 7-rays from ensemble of sources lo- 
cated at different cosmological distances. Otherwise, the 
conclusions based merely on 7-ray observations from a 
single source would be essentially model-dependent, and 



ulations concerning the spectral shape of CIB in order to 
prevent the "unreasonable" 7-ray source spectrum. This 
implies that we should accept the existence of a sharp 
pile-up in the spectrum of TeV radiation. Motivated by 
such a non-standard spectral shape of TcV-radiation, re- 
cently several extreme assumptions have been made in 
order to overcome the "IR background - TeV gamma-ray 
crisis" (Protheroe & Meyer 200C). In particular, Harwit 



therefcre would permit different interpretations concern- 



ing both the intrinsic 7-ray spectrum, and the flux of CIB. 

Presently we do face such an ambiguity, when trying 
to interpret the multi-TeV 7-ray emission of Mkn 501 
obser ved during its remarkab ly strong and long flare in 
1997 ( |Aharonian ct al. 1999b | ) . The only definite condu- 



ct al. (1999) suggested an interesting hypothesis that the 
HEGRA highest energy events are due to Bose-Einstein 
condensations interacting with the air atmosphere, and 
proposed a test to inspect this hypothesis by searching for 
peculiar features of showers detected by HEGRA in the 
direction of Mkn 501. Sub sequently, the HEGRA collabo- 



sion which can be drawn from these observations is that 
we see, most probably, significantly absorbed TeV radia- 
tion, especially at energies above 10 TeV, for which the 



ration has demonstrated ( Aharonian ct al. 200C ) that the 



optical depth could be as large as 10 ( Coppi fc Aharonian 



detected shower characteristics are in fact in good agree- 
ment with the predictions for the events initiated by or- 
dinary 7-rays. Another, even more dramatic hypothesis - 
violation of the Lorentz invariance - has been proposed by 
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analysis of the intergalactic absorption at sub- TeV and 
multi-TeV parts of the spectrum of Mkn 501 leads to two 
essentially different conclusions. The absorption-corrected 
7-ray spectrum at low energies, based on the CIB fluxes 



theorel 



reported at 2.2 and 3.5 /im (Dwek & Arendt 1998 ; Gorjian 
et al. ^000; Wright fc Johnson 200 1[) and on the current 



lem. We may add to the list of "exotic" solutions of the 
"IR background - TeV gamma-ray crisis" a less dramatic, 
in our view, hypothesis, assuming that Mkn 501 is located 
at a distance significantly less than 100 Mpc - good news 
for the advocates of non-cosmological origin of some of 
AGN and quasars (see e.g. Hoyle & Burbidge 1996; Arp 



tical predictions for NIR (see e.g. Primack et al 
200l|)7 is in a general agreement with the Synchrotron- 
self-Compton (SSC) model of X- and TeV radiation of 
Mkn 501 (see e.g. |Guy ct al . 2000; Krawczynski et al 



et al. 1997) 



2000 |; primack et al. 2001|) . On the other hand, the cor 



rections to the 7-ray spectrum at energies above 10 TeV 
based on the unexpectedly l arge CIB fluxes de t ected by 



COBE at 1 40 and 240 /im flHauscr ct al. 1998) ; |Schlegel 



et al. '. 998 ; Lagache et al. 1999 ) result in an "unreason 



able" source spectrum. Namely it implies a source spec- 
trum which sharply curves up above 10 TeV, unless we 
assume that the CIB flux at MIR between 10 and 50 
jum is quite low (a few nW/m 2 sr), and at longer wave- 
lengths it increases rapidly (yF v oc X s with s > 2) in 



order to match the COBE points (Aharonian et al. 1999b 
Coppi & Aharonian 1999b; Renault et al. 2001). Even so. 



the large DEBRA fluxes at 140 and 250 /im imply very 
flat "reconstructed" source 7-ray spectral energy distri- 
butions (SED), vS v = E 2 dN/dE = E 2 - a with a photon 
index a < 2. Such flat source spectra extending beyond 10 
TeV require, within SSC models, rather unconventional jet 
parameters, namely very large Doppler factors and very 
small magnetic fields (H. Krawczynski and P. Coppi, pri- 
vate communication; J. Kirk, private communication). 

A real trouble arises, however, when we take into ac- 
count th e recent claims about detection of CIB fl ux at 
100 /im ( Lagache ct al. 1999|; Finkbcincr ct al. 2000| ), and 
especially at 60 /im ( Finkbcincr ct al. 200C ). If we refer 
the reported fluxes to the truly diffuse extragalactic back- 
ground radiation, then a little room would be left for spec- 



Although very fascinating, it seems to us too prema- 
ture to invoke such dramatic revisions of essentials of 
modern physics and astrophysics. The nature of the FIR 
isotropic emission detected by COBE is not yet firmly es- 
tablished, and it is quite possible that the bulk of the 
reported flux, especially below 100 /im, is a result of su- 
perposition of different local backgrounds. Needless to say, 
t hat this would be the simplest solution of the problem 
( Btecker 200C ; Renault et al. 2001 ) . In this paper we adopt , 
however, that the reported FIR fluxes have universal (ex- 
tragalactic) origin. This implies that we adopt the exis- 
tence of a pronounced pile-up in the 7-ray spectrum of 
Mkn 501 above 10 TeV, but try to find an explanation of 
this spectral feature within the framework of new but yet 
conventional astrophysical scenarios. In this paper we pro- 
pose and study two potential ways to overcome the "IR 
background - TeV gamma-ray crisis" : 

(1) TeV 7-rays from Mkn 501 are not direct repre- 
sentatives of primary radiation of the source, but have 
a secondary origin, i.e. they are formed during the de- 
velopment of high energy electron-photon cascades in the 
intergalactic medium initiated by interactions of primary 
7-rays with diffuse extragalactic photons. This hypothe- 
sis implies a rather extreme assumption concerning the 
strength of the intergalactic magnetic field on > 1 Mpc, 
B < 10~ 18 G. 

(2) The "reconstructed" source spectrum of Mkn 501 
with a flat (almost constant) SED below 10 TeV, and 
a pile-up beyond 10 TeV, is partly or entirely produced 
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very limited. The only available measurement at 6 and 
15 /im on Fig. |l| derived from the ISOCAM source counts 
( Franceschini ct al. 200l| ) shown in Fig. [j] should be treated 
as a lower limit. In Fig. [I] we show also a slightly higher 
flux estimate at 15 ^m reported by |Biviano ct al. (2000 ). 
Therefore the flux estimate at the level of ~ 2 — 3nW / m 2 sr 
( Franceschini et al. 2001) as well as the lower limits based 



on the IRAS counts at 25-100 /im (Hacking fc Soifer 1991) 



Fig. 1. Cosmic background radia tion. The reporte d fluxes 
are sho wn with filled symbols: Bernstein (1999) - dia- 
monds, Wright fc Johnson (2001) - circles, IFinkbeiner 



et al. |(2000| ) - squares, [Hauser et al. (1998| ) - trian 



gles. The low limits are shown by open symbols: Pozzetti 
et a l. (1998[) - diamonds, [Biviano et al. (2000[) - trian 



do not allow firm conclusions about the depth of the MIR 
"valley" dominated by radiation of the warm dust com- 
ponent. Consequently, it does not provide sufficient infor- 
mation for definite predictions regarding the slope of the 
spectrum in the most crucial (from the point of view of 
absorption of > 10 TcV 7-rays) MIR-to-FIR transition re- 
gion. 

As long as the available measurements of CIB do not 
allow a quantitative study of the effect of absorption of 
7-rays in the intergalactic medium, we can rely only on 
model predictions or on the "best guess" shape of the CIB 
spectrum. In this regard we notice that the reported high 
FIR fluxes present a common problem for all current CIB 
models. Therefore, if one adopts that the reported FIR 
fluxes have truly diffuse extragalactic origin, an essential 
revision of the CIB models is needed in order to match the 



gle, Franceschini ct al, (2001 ) - circles, Hacking fc Soifer 
(199l|)J- squares. The CIB models are shown by solid line 



data. Such an attempt has been made recently by Primack 



et al. (2001) who showed that their semi-analytical ap- 



- Model I, by dotted line - Model II, by dashed line 
Model III, by dot-dashed line - Model IV (for details see 
the text) 



by monoenergetic ultrarelativistic beam of electrons due 
to the inverse Compton scattering in deep Klein-Nishina 
regime. In order to avoid significant radiative (synchrotron 
and Compton) losses, which otherwise would result in 
an equilibrium, E~ 2 type differential spectrum of elec- 
trons, we assume that such a beam of electrons in fact 
is a cold, conical kinetic-energy dominated ultrarelativis- 
tic wind with the bulk Lorentz factor Fo ~ 4 x 10 7 
formed beyond the accretion disk of the central black hole. 
Apparently this hypothesis implies non-acceleration origin 
of highest energy 7-rays detected from Mkn 501. 

2. Absorption of gamma-rays in CIB 

The reported fluxes and flux upper/low limits of CIB from 
optical/UV to far IR wavelengths are shown in Fig. |l|. The 
reliability and the implications of these measurements are 



proach, with a reasonable adjustment of some model pa- 
rameters, and using the empirical dust emission templates 
of Dwck ct al. (199S), can match the reported FIR fluxes. 

In Fig. [lj we show several model spectra of CIB. Since 
we are interested, first of all, in the cosmic background 
fluxes at MIR and FIR, at shorter wavelengths we adopt a 
common for all models approximation which matches the 
reported optical and NIR fluxes. The smooth template 
of the CIB spectrum in the most principal MIR-to-FIR 
transition region shown by dashed line (hereafter Model 
III) in Fig. [lj fits the reported fluxes including, within 
2a uncertainty, the 60 /jm point. Actually this idealized 
template is quite similar to the CIB spectrum shown in 
Fig. 1 of Franceschini et al. (2001 ) based on their reference 
model for IR galaxy evoluti on, as well as to som e of the 
recent theoretical models of Primack et al. (2001 ). It has 



a rather flat shape in the MIR-to-FIR transition region, 
vF u oc A s with s < 1. This results in short mean free 
paths of 7-rays above 10 TeV {L < 50Mpc; see Fig. ||), 
and consequently in a pile-up in the reconstructed 7-ray 
source spectrum (see Fig. ||) defined as: 



ME) = J ohs (E) exp [t(E)} 



(1) 



discussed in the recent review article by Hauser & Dwek 
|(200l| ). The level of the spectral energy distribution (SED) 
of CIB at optical/NIR wavelengths with the "best guess 
estimate" between 20 and 50 nW/m 2 sr is comparable with 
t he overall energy flux o f FIR of about 40-50 nW/m 2 sr 
( Pozzetti fc Madau 2000 ). This indicates that an essential 
part of the energy radiated by stars is absorbed and re- 
emitted by dust in a form of thermal sub-mm emission. 
Currently the information at mid-infrared wavelengths is 



where Jobs is the observed 7-ray spectrum, t{E) = d/ L(E) 
is the intergalactic optical depth, and d is the distance to 
the source. 

If we adopt that the reported FIR fluxes correctly de- 
scribe the level of truly diffuse background radiation, only 
two ways are left for reduction of the effect of attenu- 
ation of > 10 TeV 7-rays: (i) an ad hoc assumption of 
the CIB flux at wavelengths between 10 and 60 /im at 
the marginally acceptable (i.e. the ISOCAM low-limit) 
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Fig. 2. Mean free path of 7-rays for 4 different models 
of the CIB spectrum: solid line - Model I, dotted line - 
Model II, dashed line - Model III, dot-dashed line - Model 
IV. The horizontal lines indicate the distances to Cen A, 
M87, Mkn 501, and 3C 273 (H = 60 km/s Mpc). 



level, h ut with very rapid rise beyond 60 pm in order tc 



match ' the r 1 fluxes at FIR, and (ii) adopting for 

the Hubble constant Ho ~ 100 km/s Mpc, i.e. assuming 
the smallest possible distance to Mkn 501 d = cz/Hq = 
102 Mpc = 0.034). 

In Fig. we show 2 other model spectra of CIB (solid 
line - Model I, dotted line - Model II) which fit the data 
at NIR and FIR, but at the same time allow minimum 
intergalactic 7-ray absorption at E > 10 TeV, because 
both spectra are forced to be at the lowest possible level 
at MIR-to-FIR transition region set by the ISOCAM lower 
limit at 15 /im. 

The CIB fluxes above 140 um are appro ximated by 
the function suggested by Fixscn ct al. (1998 ). Below 140 
/xm for both I and II Models we assume the same spec- 
tral shape described by Planckian distribution. Such spec- 
tra should be considered as a lower limit, because in the 
most critical MIR-FIR transition region near 60/im it is es- 
sentially the Wicn part of the blackbody spectrum - the 
steepest (physically justified) continuous spectrum. The 

Model TT marginally qgrppg with thp 100 jrnn pnint hut nn- 



g 10 3 



Markarian 501 
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0.5 1 2 5 10 20 
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Fig. 3. Spectral Energy Distribution (SED) of Mkn 501. 
The experimental points (filled circles) correspond to the 
time-averaged spectrum of Mkn 501 during the flare in 
1997 (Aharonian ct al. 1999b). The heavy line corresponds 
to the fit to this data in the form of Eq. (0). The star 
correspond to the flux in the highest energy bin around 
21 TeV obtained after the reanalysis of the same data 



set but with improved energy resolution (Aharonian et al 



2001b). The vertical line at 17 TeV indicates the edge 



of the spectrum of Mkn 501 measured by HEGRA with 
high statistical significance. The solid, dotted, dashed and 
dot-dashed lines represent the reconstructed (absorption- 
corrected) spectra of 7-rays for the CIB Models I, II, III 
and IV, respectively (Hq — 60 km/s Mpc). 



an absolute value for M kn 501 t(E) ps 0.16(£7l TeV) 
1999b|) . Therefore, in the absorption- 



( [Aharonian et al. 

corrected (source) spectrum of 7-rays given by Eq. 
the optical depth r almost compensates the exponential 
term of the observed spectrum of Mkn 50lQ, which from 
0.5 TeV to 20 TeV can be prese nted in a simple "power-l aw 
with exponential cutoff" form ( Aharonian et al. 1999b| ) 



Jobs(E) cx E exp (-£76.2 TeV) 



(2) 



Recently, the HEGRA collaboration publis hed the re- 
sults of reanalysis of the spectrum of Mkn 501 (Aharonian 



rWpgtll-natog thp flnv at F.CI jim hy a faptnr nf K rnmparprl 

with the flux reported by Finkbeiner et al. (2000). This 
model significantly suppresses the photon density between 
50 and 100 fim, and correspondingly allows larger mean 
free paths for 7-rays with energy more than several TeV 
(Fig. |). Such a SED of CIB results in almost E~ 2 type 
"reconstructed" spectrum of 7-rays from Mkn 501 at en- 
ergies between approximately 2 and 20 TeV (Fig. ||). This 
has a simple explanation. For the constant SED of CIB 
(the flat part of the spectrum 2 in Fig. |l|), ncBPi(e) oc e~ 2 , 
the photon-photon optical depth is proportional to E with 



ct al. 2001b) using an improved method of shower energy 



reconstruction with 10 to 12 per cent energy resolution 
(Hofmann 2000). Apparently the improved energy resolu- 
tion is more relevant for derivation of the shapes of steep 
7-ray spectra, e.g. spectra with exponential or sharper cut- 
offs. The new HEGRA spectral analysis confirmed the re- 
sults of the previous study up to energy 17 TeV, but shows 



Note that somewhat different CIB models with more 
complex spectral shapes can also lead to flat, E~ 2 type 



power-law absorption-corrected 7-ray spec t ra (see e.g. Coppi 
^ Aharonian 1999b| ; [Konopelko et al. 199S| ; |Guy et al. 2000] ) 
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a steeper spectrum beyond that energy ( Aharonian ct al 
2001b| Therefore in Fig. we show only one point of the 
new analysis - the flux in the highest energy bin around 21 
TeV (star) which lies significantly below the flux estimate 
of the previous study (open circle). We should notice, how- 
ever, that the points above 17 TeV in both papers do not 
have high statistical significance anyway, therefore below 
we limit our study only by the energy region up to 17 
TeV (i.e. by the region on the left side of the vertical line 
shown in Fig. |3|). In this energy interval the energy spectra 
obtained by new and old methods are perfectly described 
(|Aharonian ct al. 2001b] ) by Eq. (§). 

Apparently, for conventional 7-ray production mecha- 
nisms the E~ 2 type 7-ray source spectrum sounds more 
comfortable than the spectrum containing a sharp pile- 
up. Therefore, the Model II in Fig. [I] can be treated as the 
"most favorable" CIB spectrum, although its flat shape 
in the broad MIR-to-FIR transition region with very fast 
increase beyond 60 /im hardly match the current theo- 
retical descriptions of CIB, as well as the FIR spectra 
observed from nearby galaxies. Such a comparison is of 
special interest, especially because the bulk of the inte- 
grated extragalactic background radiation has been gen- 
erated presumably at z < 1 ( Harwit 1999| ; Franccschini 



et al. 2001). In particular, the shape of the Model II con- 
tradicts to the average SED of the ISOCAM sources which 
not only contribute a dominant fraction to CIB at MIR, 
but likely are major contributors at longer wavelengths as 
well ( Franceschini et al. 2001). 

All current CIB models have a problem to accommo- 
date the reported flux at 60 fim. If this point, however, 
is the representative of the truly diffuse flux, we must as- 
sume a CIB spectrum close to the Model - 1 shown in Fig. [j] 
by the solid line. This Model assumes a very steep slope 
between 30 and 60 /im. A steeper spectrum (i.e. steeper 
than the Wien tail of the black-body radiation) in this nar- 
row wavelength band hardly could be physically justified. 
Moreover, already 60 /xm photons have sufficient energy 
for effective interaction with > 10 TeV 7-rays, therefore 
we cannot suppress anymore the severe 7-ray absorption 
by speculating about the spectral shape of CIB (see Fig. |] 
and the related discussion below). 

Finally, in Fig. [j] we show one more possible model 
spectrum of CIB (dot-dashed curve - Model IV) which as- 
sumes significantly, by a factor of 2.5, higher flux at 15 fim 
compared to the reported ISOCAM lower, and smoothly 
passes through the low edges of the error bars of reported 
fluxes at 100 and 140 /im. Surprisingly such a high at MIR 
spectrum does not result in an unusual 7-ray spectrum 
as long as it concerns the energy region below 10 TeV. 
Namely, at these energies we obtain an almost single hard 
power-law 7-ray source spectrum with photon index less 
than 2. Above 10 TeV we again observe a pile- up which 
however in this case is less pronounced than in the case of 
Models I and III. 

A 7-ray photon with energy E propagating trough 
isotropic photon field can interact, via electron-positron 
pair production, with ambient photons of energy e > 




A (/Am) 

Fig. 4. The contour map of the function k(X,E). The 
solid, dotted and dashed lines show the levels of n(\, E) 
for the CIB Models I, II, and III, respectively. The heavy 
solid line represents the threshold of 77 pair production. 



e th = (m c c 2 ) 2 /E ~ 0.26(-B 7 /lTeV)- 1 eV or of wavelength 
A = 4.8C£f 7 /lTeV) /im. The latter relation is represented 
in Fig. [| by the solid line. The cross-section of 77 interac- 
tions averaged over the directions of background photons 
peaks at e max ~ 3.5e t hrcsh with a™ x ~ (l/4)cr T (see, e.g. 
IVassiliev 2000|) 

Therefore, even for a broad, e.g. power-law spectrum 
of background radiation, the most contribution to the op- 
tical depth r comes from a narrow spectral band of CIB 
within e max ± Ae with Ae ~ l/2e max . For typical spectra 
of CIB with two distinct NIR and FIR bumps and a MIR 
valley, the "7-ray energy-CIB photon wavelength" relation 
is more complicated, namely the relative contributions of 
different spectral intervals of CIB to the optical depth To 
significantly depend on the energy of the primary 7-ray 
photon. It is convenient to describe this dependence by 
the ratio k(E, A) = t(E, \)/tq(E), where 



t(E,X) 



1 (E 1 , A)ncBR(A) dA , 



(3) 



and To = t(E~ ( ,X = 0) is the optical depth integrated 
over the entire spectral range of CIB nciB(A) above the 
threshold A th ~ 4.8 (J3 7 /l TeV) /mi. 

In Fig. |l| the contour map of the function n{E, A) is 
shown for 3 different models of CIB in the spectral range 
of 7-rays from 1 to 100 TeV, and for the CIB photons 
from 1 to 500 /zm. For the given energy of 7-ray photon 
Eq , the spectral region of CIB responsible for the fraction 
£ of the total optical depth To is [A^,Ath], where A^ and 
Ath are the abscissas of the points where the horizontal 
line E = Eq intersects the corresponding level curve, and 
the threshold line, respectively. In Fig. |J four levels for 
£ = 10,50,90, and 99%, are shown. By definition, £ = 
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corresponds to the threshold boundary £?TeV = 0.21 A p . 
It is seen that at Eq = 17 TeV approximately 50 per cent 
of the total optical depth tq is contributed by background 
photons with wavelengths longer than 50/im. Therefore 
even the sharp cutoff of the CIB flux below 50/im cannot 
prevent the large optical depth of 17 TeV 7-rays, as far as 
we accept that the 60/im flux has a cosmological origin, 
and reflects the level of CIB at these wavelengths. This 
is the case of Model I for which the free path of 17 TeV 
7-rays is about 20 Mpc, and therefore for the distance to 
the source of 170 Mpc (Ho = 60 km/s Mpc) the optical 
depth To ~ 8.5. The corresponding absorption factor is so 
large (exp(— 8.5) ~ 2 x 10~ 4 ) that the sharp file-up in the 
absorption-corrected spectrum becomes unavoidable. 

Since the exponent in the absorption factor is very 
large, even a relatively small reduction of the optical depth 
To may help to reduce significantly the effect of 7-ray at- 
tenuation. In particular, the Model III, which has a more 
realistic spectral shape at mid infrared wavelengths, and 
fits the FIR data above 100 /im, allows a bit larger (25 per 
cent or so) free paths at highest energy 7-rays. Although 
this makes the pile-up less pronounced, but still does not 
eliminate it altogether. Only the CIB Model II gives a 
comfortable, E~ 2 type spectrum. 

We may reduce the absorption effect furthermore, as- 
suming a larger, although presently less favored value for 
the Hubble constant, Ho = 100 km/s Mpc. The latter 
makes smoother the pile-up for the Model I, and almost 
removes it for the Models III and IV. Therefore we may 
conclude that it is possible, in principle, to avoid the "dis- 
turbing" sharp turn-up in the intrinsic 7-ray spectrum of 
Mkn 501, if we adopt for CIB a model like the Model III, 
and assume a very large value for the Hubble constant. 

Below we discuss, however, possible solutions which 
allow accommodation of both higher CIB fluxes and more 
realistic value for the Hubble constant. 

3. The effect of cascading in the intergalactic 
photon fields 

Generally, the propagation of high energy 7-rays trough 
a low frequency photon field cannot be reduced to the 
.effect of 77 absorption, When a 7-ray is absorbed 



hundred TeV interact mostly with the infrared and opti- 
cal background photons. The inverse Compton scattering 
of electrons does not have kinematic threshold, therefore 
the electrons interact predominantly with much denser 
CMBR. 

For intrinsic 7-ray spectra harder than E~ 2 extend- 
ing to energies E > 100 TeV , the cascade spectrum may 
well dominate over the primary 7-ray spectrum. On the 
other hand, the cascade spectrum typically has a stan- 
dard shape which slightly depends on the primary 7-ray 
spectrum. This makes rather attractive the idea of inter- 
pretation of one of the most remarkable features of the 
TeV radiation of Mkn 501, namely the surprisingly stable 
spectral shape of the source in high state, despite dramatic 
variation of the absolute flux on timescales less than sev- 
eral hours (Aharonian ct al. 1999b a]). Moreover, as it was 
noticed by 



Aharonian et al. 



(1999bQ , 



the shape of the 

TeV spectrum of Mkn 501 given by Eq. (||) resembles the 
7-ray spectrum formed during the cascade development in 
the photon field. Therefore, the quantitative study of this 
effect presents a definite interest. 

The interpretation of the observed TeV emission from 
Mkn 501 in terms of the intergalactic cascade radiation 
requires an extremely low intergalactic magnetic field. 
Indeed, for the field exceeding 10~ 12 G, the 7-rays of the 
cascade origin should be observed in a form of an extended 
emission from a giant isotropic pair halo with an angular 
radius more than several degree (Aharonian et al. 1994). 
Both the detected angular size and the time variability of 
the TeV radiation from Mkn 501 excludes such a possi- 



bility (Aharonian ct al. 1999a, 2001c). For lower magnetic 
fields, the cascade 7-rays penetrate almost on a straight 
line, thus the first argument based on the detected angu- 
lar size becomes less stringent. Note however that for the 
distance to the source d > 100 Mpc even tiny deflections 
of the secondary (cascade) electrons by the intergalactic 
magnetic field would lead to non-negligible t ime delays of 
arriving 7-rays (|Plaga 1995| ; |Kronbcrg 1995|) : 



At B ^ 10(d/170Mpc)(£yiTeV)" 2 (5/10- i8 Gr h 



(4) 



This would smear out the time variation of TeV fluxes 



simple 



its energy in fact is not Inst The secondary electrons and 



positrons create new 7-rays via inverse Compton scatter 
ing; the second generation 7-rays produce new (e + ,e _ ) 
pairs, thus an electromagnetic cascade develops. Actually, 
in the intergalactic space this process is inevitable, and 
it may significantly contribute to the isotropic (extra- 
galactic) 7-ray background radiation (Protheroe & Stanev 
1996; Coppi & Aharonian 1997). Note that in the energy 



observed on timescales less than several hours ( Aharonian 
ct al. 19994 , unless B < 10~ 18 G§ 

For such a low intergalactic magnetic field, we have 
studied the spectral properties of electromagnetic cas- 
cades generated by the primary multi-TeV radiation of 
Mkn 501 in the intergalactic medium, using a fast numer- 
ical method based on the solution of adjoint cascade equa- 



region of interest the 7-rays and electrons interact with 
photon fields of different origin. Although the energy den- 
sity of the 2.7 K cosmic microwave background radiation 
(CMBR) well exceeds the density of other ("starlight" and 
"dust") components of the cosmic background radiation 
(see Fig. [j]), because of the kinematic threshold of the ra- 



tions (Uchaikin & Ryzhov 1988). Two processes have been 
taken into account: electron-positron pair production and 
inverse Compton scattering. The details of calculations 
will be published elsewhere. 



Actually, due to non-zero (r 



action 77 



e the 7-rays with energy less than several 



,c 2 /E) emission angles of 
the secondary products in the 77 — » e + e~ and ej — > ej' reac- 
tions, we should expect a non-negligible broadening, and con- 
sequently time delays of the cascade radiation even at the ab- 
sence of magnetic field ( Sheng & Cheng 199C ) . 
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Fig. 5. The cascade development initiated by primary mo- 
noenergetic 7-rays with energy (a) Eq = 10 4 TeV and (b) 
E = 10 2 TeV for the CIB Model I. The curves show the 
number of cascade 7-rays in different energy intervals as 
a function of the penetration depth. 
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Fig. 6. The same as in Fig. S but for the power-law 
spectrum of primary 7-rays, Jq(E) oc E~ 2 extending to 
E max = 10 2 TeV (dashed lines) and E max = 10 4 TeV 
(dot-dashed lines). The solid lines show the evolution of 
the number of un-absorbed primary 7-rays, i.e. without 
the contribution of the secondary (cascade) component. 



The idea behind the attempt to interpret TeV radia- 
tion of Mkn 501 by the intergalactic cascade is the follow- 
ing. The number of photons in a given energy band of the 
cascade spectrum may essentially exceed the number of 7- 
rays in the absorbed primary 7-ray spectrum. The 7-rays 
of highest energies, during their propagation trough the 
intergalactic photon fields, initiate cascades, which pro- 
duce many lower energy photons, i.e. transfer the energy 
from the primary highest energy photons to lower energy 
7-rays. Thus, the cascade somewhat "move" the source 
closer to the observer. In order to illustrate this effect, in 
Figs. ^ and ^ we show the number of cascade 7-rays in 
5 energy bands as a function of the penetration depth. In 
the case of a cascade initiated by a monoenergetic primary 
7-ray of energy Eq (Fig. [^), the number of photons in all 
energy bands sharply increases and reaches its maximum 
at a distance i?* which is determined approximately by 
the condition t(E ) ~ 1 or i?„ ~ L(E ), where L(E ) is 
the mean free path of primary photons. Beyond -R* the 
cascade develops more slowly. This stage is characterized 
by a competition between production and absorption pro- 
cesses. While the low energy (0.1 - 0.3 TeV) 7-rays con- 
tinue to (slowly) grow, the number of higher energy (10-30 
TeV) 7-rays drops beyond 10 Mpc, the reason being the 
expiration of the "fuel" , i.e. particles of sufficient energy, 
which could support further development of the cascade 
at these energies. It is interesting to note that the overall 
picture slightly depends on energy of the primary photon 
(compare curves in Fig. || calculated for E — 10 4 and 
10 2 TeV). This explains why the spectrum of 7-rays of 



the well-developed cascade becomes almost independent 
of the primary 7-ray spectrum. 

The case of the cascade initiated by a broad-band spec- 
trum of primary 7-ray photons is more complicated be- 
cause the observer detects a mixture of the primary (un- 
absorbed) and secondary (cascade) 7-rays. Two example 
of cascades triggered by a E~ 2 type primary 7-ray spec- 
trum with £ max = 10 2 TeV (dashed lines) and E max = 10 4 
TeV (dot-dashed lines) are shown in Fig. ^. In order to 
demonstrate the photon excess caused by the cascade de- 
velopment, we show also the evolution of the number of 
photons calculated for the simple absorption effect, i.e. 
J(E) = J exp[-R/L(E)]. 

In Fig. we show the spectra of cascade 7-rays 
from Mkn 501, calculated for the Hubble constant Hq = 
60 km/s/Mpc. The spectrum of the primary 7-radiation 
was assumed to be power-law Jq(E) cx E~ a with a = 1.9. 
For softer source spectra (a > 2), the cascade has a small 
impact on the resulting spectrum. For harder primary 
spectra with a < 2, the most power is radiated in the 
high energy range, and therefore the resulting spectrum is 
dominated by the "cascade" component. 

For a narrow, e.g. Planckian type spectrum of tar- 
get photons with temperature T r , the spectrum of well- 
developed cascade 7-rays has a standard power-law shape 
with photon index 1.5 at energies E <C E* ~ m 2 c 4 /fcT r , 
and sharp cutoff beyond E* . For a broad-band background 
photon field, the spectral shape of CIB plays an essential 
role in forming the cascade 7-ray spectrum, especially in 
the region of the cutoff. This can be seen in Fig. [?] where 



8 



F.A. Aharonian et al.: Highest energy gamma-rays from Mkn 501 



1 n2 
10 


I 


i 






w 






w 






'a 

o 






CD 






CM 






o 10 1 






m 






V 

1 \\ 


> 


\ 






' V 






Markarian 501 








1 


i ■ ■ 










0.5 1 3 5 10 30 

E (TeV) 



Fig. 7. Intcrgalactic cascade radiation spectra initiated 
by primary 7-rays from Mkn 501. The solid, dotted, and 
dashed lines correspond to the Model I, Model II, and 
Model III of CIB, respectively. The light and heavy lines 
represent the predicted cascade 7-ray fluxes before and 
after convolution with the Gaussian type energy spread 
function with 20% energy resolution. 
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Fig. 8. Intergalactic cascade radiation spectra initiated by 
primary 7-rays from Mkn 501. The solid, dashed, dot- 
dashed and dotted curves are the cascade 7-ray spectra 
corresponding to the the photon indices of the primary 7- 
radiation a— 1.5, 1.8, 1.9 and 2, respectively. The curves 
are obtained after the convolution of the cascade 7-ray 
spectra with the Gaussian type energy spread function 
with 20% energy resolution. All curves are calculated for 
the CIB Model III. 



the 7-ray spectra are calculated for 3 different models of 
CIB presented in Fig. |[ For comparison, we show also the 
fluxes of Mkn 501 as measured by HEGRA during the high 
state of the source in 1997. The experimental points (filled 
circles up to 17 TeV and an open circle at ~ 21 TeV) are 
obtained with an energy resolution of about 20%. In addi- 
tion, they contain w 15% uncertainty on the energy scale 
throughout the entire interval from 0.5 to 20 TeV, and up 
to factor of 2 statistical and systematic uncertainties in 



flux estimates above 10 TeV (Aharonian et al. 1999b) 

Within these uncertainties, the cascade 7-ray spec- 
trum corresponding to the CIB Model II (light dotted line) 
agrees with the HEGRA points, while the 7-ray spectra 
calculated for the CIB Models I and III (light solid and 
dashed curves, respectively) pass significantly below the 
measured fluxes at E > 10 TeV However, at these energies 
the spectra are very steep, therefore we must to take into 
account the limited instrumental energy resolution when 
comparing the theoretical predictions with the differen- 
tial flux measurements. Namely, the experimental fluxes 
should be compared with the predictions after convolving 
the theoretical 7-ray spectra with the energy resolution 
function. In Fig. [?] the corresponding curves, obtained as- 
suming Gaussian type instrumental energy spread func- 
tion G{E - E') = exp(-(£ - E') 2 /2a 2 )/{V2^o-) with 
constant a — E/5, i.e. 20% energy resolution, are shown 
by heavy solid (Model I), dotted (Model II) and dashed 
(Model III) lines. It is seen that after this procedure the 
7-ray spectra corresponding to the CIB Models I and III 



quite satisfactorily fit the HEGRA points. At the same 
time the CIB Model II gives noticeably higher 7-ray fluxes. 

The HEGRA stereoscopic system of telescopes has a 
potential to reduce the energy resolution down to ~ 10% 
(Hofmann 2000). The recent re-analysis of the 1997 high 
state Mkn 501 data (Aharonian et al. 2001b) based on 
the improved energy reconstruction method confirmed the 
previous flux measurements, except the flux upper limit 
at ~ 21 TeV which appeared by a factor of 2.5 below the 
previous flux estimate (the star and open circle in Fig. ^, 
respectively). Therefore the 20% energy resolution used 
for calculations presented in Fig. ^should be considered as 
a rather conservative assumption. Nevertheless this seems 
to be, an appropriate assumption, at least for the purposes 
of this paper, given the lack of independent experimental 
methods for cross-calibration of the energy resolution of 
ground-based Cherenkov telescopes. 

At high energies the spectrum of cascade photons 
weakly depends on the the power-law index a and the 
maximum energy i? max of primary 7-rays, provided that 
a < 2 and E max ^ 10 2 TeV. In Fig. || four cascade 7-ray 
spectra are shown calculated for a = 1.5,1.8,1.9,2 and 
£ max = 10 3 TeV and assuming 20% detector energy res- 
olution. All spectra are calculated for the CIB model III. 
The photon indices of primary 7-rays a — 1,8, 1.9, and 2 
provide very similar fits for the HEGRA points through- 
out the entire energy interval from 0.5 TeV to 20 TeV. 
They become distinguishable at low energies, but mainly 
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because of the contributions of primary (un-absorbed) 7- 
rays which dominate below 1 TeV. The very hard primary 
spectrum with a — 1.5 results in a flatter cascade 7-ray 
spectrum compared to the measurements. On the other 
hand, the primary photon index a should not exceed 2.0- 
2.1, otherwise the efficiency of the cascade contribution be- 
comes very small. Hoverer, for the photon index of about 
2 the primary spectrum has to extend beyond 100 TeV in 
order to provide an adequate contribution from the cas- 
cade component. 

4. Gamma-radiation produced by ultrarelativistic 
unshocked jets 

The " intergalactic cascade" hypothesis requires extremely 
small intergalactic magnetic field of about B ~ 10~ 18 G 
and extension of the primary 7-ray spectrum well beyond 
100 TeV. In the case of failure of any of these two condi- 
tions, the contribution of the secondary 7-rays produced 
in the intergalactic medium becomes negligible, and thus 
the attenuation of 7-rays is reduced to the simple absorp- 
tion effect. Consequently, the absorption-corrected spec- 
trum of 7-rays given by Eq. (Q) adequately represent the 
intrinsic TeV spectrum of the source. If so, we may face a 
difficulty with the current common belief that the source 
spectrum of 7-rays should have a "decent" shape, i.e. be 
in accord with the predictions of conventional astrophys- 
ical scenarios suggested for blazars. Indeed, none of the 
models of blazars, in general, and of Mkn 501 in particu- 
lar, allow the striking feature which appears unavoidably 
in the absorption-corrected spectrum of Mkn 501 above 
10 TeV (see Fig. ||), if we adopt that the reported high 
FIR fluxes are due to the truly extragalactic background 
radiation. All versions of both the electronic and hadronic 
models, suggested so far, predict smooth broad-band 7-ray 
spectra with characteristic (gradual or sharp) steepening 
above 10 TeV. 

Formally, it is possible to reproduce a 7-ray spectrum 
with a sharp pile-up assuming very narrow features in the 
distribution of accelerated particles. It should be noticed 
in this regard that almost all particle acceleration models 
predict power-law distributions with high energy cut-offs. 
While the cutoff energy Eq can be estimated, generally 
quite confidently, from the balance between the particle 
acceleration and the energy loss rates, the shape of the re- 
sulting particle spectrum in the cutoff region depends on 
many specific acceleration and energy dissipation mecha- 
nisms. Remarkably, even within the "ordinary" shock ac- 
celeration scenarios we may expect not only spectral cut- 
offs, but also pronounced pile-ups preceding the cutoffs 
(e.g. [Malkov (19971); [Melrose fc Crouch (1997Q ; |Protheroe 
fc Stanev (1999Q ; Drury et al. (1999Q , see however [Kirk 




(2001)). Within the proton-synchrotron model of blazars, 



such a pile-up would result in a b ump in the sync hrotron 
TeV emission at ~ 0.3r}- 1 D i TeV (|Aharonian 2000| ), where 
77 > 1 is the so-called gyrofactor. Thus, the location of the 
bump at E > 17 TeV should require very large Doppler 
factor Dj > 50rj. In the leptonic models the pile-up in 



Fig. 9. Inverse Compton spectrum of cold unshocked ul- 
trarelativistic jet with bulk motion Lorentz factor T = 
3.33 x 10 7 . The radiation component associated with 
comptonization of ambient optical photons with narrow 
spectral distribution is shown by the dashed line I. The 
IC spectrum on ambient photons with broad-band spec- 
tral distribution is shown by the dashed line II. The heavy 
solid line represents the superposition of these two com- 
ponents. Formally, the dashed line II could be treated as 
residual of the total TeV source emission after subtraction 
of the unshocked wind component I, and therefore can be 
referred to the IC radiation of blobs in shocked jet (see 
Fig. [ll|). Fits to the observed flux of Mkn 501 are shown 
by thin solid lines. Curve 1 corresponds to the fit given 
by Eq. @) and curve 2 - to the steepest possible spectrum 
above 17 TeV based on the re analysis of Mkn 501 HEGRA 
data ( Aharonian et al. 2001b ). Absorption-corrected spec- 
trum of Mkn 501 for CIB Model I is shown by dotted lines 
1' and 2' for the fits to the observed spectrum 1 and 2 re- 
spectively. 



the electron spectrum would also result in corresponding 
features in the synchrotron X-ray and inverse Compton 
TeV 7-ray spectra. The absorption-corrected spectra of 
Mkn 501 shown in Fig. [^ require very sharp pile-up in 
the electron spectrum. Whether such a pile-up could be 
formed in realistic particle acceleration models of small 
(sub-pc) jets of blazars is a question of future detailed 
studies. 

In this paper we offer a different, non- acceleration 
scenario which postulates that the TeV radiation of 
Mkn 501 is a result of comptonization of the ambient 
low-frequency radiation by ultrarelativistic jet-like out- 
flow with a Lorentz factor of the bulk motion of about 
r ~ (3 - 4) x 10 7 . 

The relativistically moving plasma outflows in forms 
of jets or winds, are common for many astrophysical phe- 
nomena on both galactic or extragalactic scales (see e.g. 
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Mirabel 2000). Independent of the origin of these rela- 
tivistic outflows, the concept of the jet seems to be the 
only successful approach to understand the complex fea- 
tures of nonthermal radiation of blazars, microquasars 
and GRBs. The Lorentz factor of such outflows could be 
extremely large. In particular, in the Crab Nebula the 
Lorentz-factor of the MHD wind is estimated between 
10 6 and 10 7 (|Rees k Gunn 1974| ; |Kundt k Krotschcck 
1980; Ifcnncl k Coroniti 1984 ). The conventional Lorentz- 
factors of jets in the inverse Compton models of 7-ray 
blazars, are rather modest, T ~ 10. However there are 
no apparent theoretical or observation arguments against 
the bulk m otion with extreme L orentz-factors (see e.g . dis- 
cussion by 



Celotti et al. (199S| ) on Mkn 421). |Meszaros k 



Rees ( .997 ) have shown that in the context of cosmolog- 
ical GRBs the magnetically dominated jet-like outflows 
from stellar mass black holes may attain extreme bulk 
Lorentz factors exceeding 10 6 . We are not aware of simi- 
lar calculations (and conclusions) concerning the massive 
black holes - the engines of AGN. Therefore we will limit 
our discussion by postulating existence of an outflow in 
Mkn 501 which moves, at least at the initial stages of its 
propagation in the vicinity of the central black hole, with 
an extreme bulk Lorentz factor > 10 7 . If true, it is almost 
obvious that the outflow should have MHD origin, the en- 
ergy being extracted from the rotating black hole through 



e.g. the Blandford-Znajek type mechanism (see Blandford 
(2000j)]for a recent review). In proximity of the accretion 
disk, the outflow should be Poynting-flux dominated in 
order to avoid the Compton drag, and only at relatively 
large distances from the central object, where the photon 
density is significantly reduced, an essential fraction of the 
electromagnetic energy must be transfered to the kinetic 
energy of bulk motion. 

From a cold relativistic outflow we do not expect sig- 
nificant synchrotron radiation. Indeed, although the en- 
ergy of electrons in the frame of observer can be as 
large as 10 TeV, they move together with magnetic 
field, and thus they do not emit synchrotron photons. 
Nevertheless, the cold ultrarelativistic outflow could be 
visible through the inverse Compton 7-radiation of wind 
electrons. Apparently, in this scenario the dense target 
photon field is the second important ingredient for effec- 
tive production of 7-rays. Remarkably, even in the case 
of relatively weak BL Lac objects there could be sev- 
eral important sources of infrared and optical emission 
within the inner sub-parsec region of the central source 
- IR emission from the dust torus, broad-line emission 
from fast moving clouds, starlight, etc. It is easy to show 
that the Compton optical depth tq in this region could 
be as large as 1 ( Celotti ct al. 1998| , see e.g.). In fact, 
we have to introduce some additional conditions to avoid 
the destruction of the jet due to the Compton drag. On 
the other hand, since for tq < 1 the 7-ray luminosity is 
proportional to the product Lj x re, the optical depth 
should be close to 1 in order to avoid assumptions which 
would require very large jet power L y Obviously, the 
most favorable value for tq < 1 lies between 0.1 and 1. 




Fig. 10. The same as Fig.g, but for CIB Model III. 



Because of extremely large bulk Lorentz factor T > 10 7 , 
the Compton scattering on the ambient optical photons 
and near-infrared photons proceeds in deep Klein-Nishina 
regime. Therefore the resulting 7-radiation should have a 
very narrow distribution with a maximum at the edge of 
the spectrum E w E e = m c c 2 r. Meanwhile, the IC scat- 
tering on low-frequency target photon components, e.g. 
on far infrared photons of the dust torus, would still take 
place in the Thomson regime, and thus would produce a 
smooth broad-band 7-ray spectrum. 

In Figs. I and [l(] we demonstrate that the overall 
absorption-corrected spectrum of Mkn 501 can be satis- 
factorily interpreted in the terms of bulk inverse Compton 
emission of the jet, assuming a specific ambient radiation 
field consisting of two - narrow (Planckian) type radiation 
with temperature kT ~ 1 eV (curve I) and broader IR ra- 
diation, which formally could be presented in a power-law 
form n(e) oc e _m (curve II). More specifically, the best 
fit for the CIB Model I is achieved for kT = 2 eV and 
m = 1.8, with the ratio of the energy densities of two 
components of about r — w±/w2 — 163. The CIB Model 
III requires similar, although slightly different parameters: 
kT = 0.2 eV, m — 1.8, and r = 9.5. In both cases the bulk 
Lorentz factor T = 3.33 x 10 7 is assumed. 

The "Klein-Nishina" component, caused by inverse 
Compton scattering on the "hot" narrow photon distri- 
bution, is of prime interest because this component pro- 
vides the most critical part of the spectrum - the pile-up. 
Moreover, we can omit the second ("Thomson") compo- 
nent of 7-radiation (for example assuming larger values 
for the ratio r = w\/w2 compared to the ratios used in 
Figs. H and |Io| ). This should not be a problem for the 
model, because we may compensate this loss by assum- 
ing that the low (E < 10 TeV) energy part of the 7-ray 
spectrum is formed at later stages of evolution of the jet, 
namely due to transformation of a fraction of the bulk 
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kinetic energy of the jet into a random energy of rela- 
tivistic electrons, e.g in the form of so-called blobs (plas- 
mons) moving with Lorentz factor of about 10. This can 
be realized through the acceleration of electrons by the 
terminal shock which may be created at interaction of the 
jet with the ambient matter (here again, we could see an 
interesting analogy with the Crab Nebula). Actually this 
is the stage which is treated in all conventional, SSC or 
external Compton models of X-ray and 7-ray emission of 



TeV blazars (Bikora fc Madejski 2001) with "standard" 



magnetic field of about 0.1-1 Gauss ( 


Inouc & Takahara 


199£; 


Jednarck & Protheroe 1997; 


Tavecchio et al. 1998; 


Mastichiadis & Kirk 1997; 


Krawczynski et al. 2000 


). The 



5. Discussion 

It is believed that the diffuse extragalactic infrared back- 
ground radiation may have a dramatic impact on the mod- 
els of TeV blazars. In particular, the recently reported high 
CIB fluxes, both at NIR and FIR bands, imply that we 
detect significantly absorbed TeV radiation even from rel- 
atively nearby objects like Mkn 421 and Mkn 501. Due 
to the energy-dependent mean free path of 7-rays in the 
intergalactic medium, the detected spectra of TeV emis- 
sion from extragalactic objects at cosmological distances 
significantly deviates from the source spectra. The inter- 
galactic absorption does not simply imply spectral cut- 
offs, but rather modulates the primary spectrum from 
sub- TeV to multi-TeV energies. For some specific CIB 
models, the "high-state" TeV spectrum of Mkn 501 after 
correction for the intergalactic absorption becomes rather 



flat (Fig. H ), dN/dE cx E~ 2 (e.g. 


Coppi & Aharonian 


1999a; 


Konopclko et al. 1999|) or even harder ( 


Guy et al. 


200C ; 


Renault et al. 2001 


; ie Jager & Stecker 2001), 



if the CIB density o f the mid-infrared "valle y" exceeds 
vF v > 4 nW/m 2 sr ( [Aharonian et al. 1999b|) . Although 
formally such spectra can be described within the conven- 
tional SSC models, the detailed treatment of the problem 
within the one-zone SSC model requires quite uncomfort- 
able jet parameters, in particular very large Doppler factor 
Sj > 100, and very weak magnetic field B < 0.01 G. If so, 
this perhaps would require a serious consideration alterna- 
tive models for TeV emission like the proton-synchrotron 
model which can satisfactorily fit the data of Mkn 501 and 
Mkn 421 even for quite high CIB fluxes (Aharonian 2000). 

However all these models fail to explain the sharp pile- 
up which appears at the end of the "reconstructed" spec- 
trum of Mkn 501, if the reported FIR fluxes correctly de- 
scribe the level of the diffuse cosmic background radiation. 
Remarkably, such a pile-up is unavoidable not only due to 
the reported extremely high flux at 60/im, which obvi- 
ously needs further confirmation, but also, albeit in a less 
distinct form, due to the fluxes at longer, A > 100 ^m 
wavelengths reported by several independent groups. To 
avoid such an unusual shape of 7-radiation, recently sev- 
eral dramatic assumptions have been proposed concerning 



the nature of the detected signals or the violation of the 
Lorenz-invariance at very high energies. 

The main objective of this work was the attempt to 
investigate other possible ways to overcome the "IR back- 
ground - TeV gamma-ray crisis" . In particular, we studied 
two possibilities: (i) assuming that the observed TeV spec- 
trum of 7-rays from Mkn 501 is formed due to the cascade 
initiated by primary 7-rays in the intergalactic medium, 
and (ii) assuming that the detected highest energy 7-rays 
are result of bulk-motion comptonization of the ambient 
optical photons by the ultrarelativistic unshocked conical 
wind emerging from the central source. 

5.1. Rectilinear intergalactic cascade radiation 

The electromagnetic cascades in the extragalactic pho- 
ton fields result in a non-negligible increase of the ef- 
fective propagation path of 7-rays. Remarkably, the re- 
sults presented in Sec. 3 show that the cascade 7-ray 
spectra can quite well fit the time-averaged high state 
spectrum of Mkn 501. An important feature of the well- 
developed cascade radiation is its standard spectral shape 
which is determined by the propagation effects and weakly 
depends on the details of the primary radiation of the 
central source. Thus, this hypothesis could give a nat- 
ural explanation for the quite stable spectral shape of 
the source in high state, despite dramatic variation of 
the absolute flux on timcscalcs less than several hours 



( Aharonian et al. 1999b , a|). Note that the spectral change 
during the stro ng April 16, 1997 flare repo rted by the CAT 
collaboration ( Djannati-Atai et al. 1999 ), as well as the 
noticeable steepening of the spectrum in a low state of 



Mkn 501, found by the HEGRA collaboration (Sambruna 
^t al. 200C ; Aharonian et al. 2001a) do not contradict to 
this statement. These effects, in fact, could be caused by 
variations of the ratio of the "cascade" component to the 
overall ( "unabsorbed" plus "cascade" ) flux of 7-rays due 
to, for example, reduction of the maximum energy in the 
primary 7-ray spectrum at the low or quiescent states. 

The "intergalactic cascade" hypothesis requires ex- 
tremely small, at least in the direction of Mkn 501, in- 
tergalactic magnetic fields, in order to avoid the signifi- 
cant time delays. Even a tiny intergalactic magnetic field 
of about 10~ 18 G should lead to delays of arrival of TeV 
7-rays, compared to the associated low energy photons, 
by more than several days. This obviously contradicts 
the observed X/TcV correlations (e.g. [Plan et al . 199*8; 



|Catanese et al. 1997 ; Krawczynski et al. 2000[) . Although 



indeed very speculative, formally such small fields can- 
not be excluded, especially if we take into account that 
the typical scale of the so-called intergalactic voids, where 
the magnetic field could be arbitrarily small, is estim ated 
as - 120(i? /100 km/s/Mpc)" 1 Mpc flEinasto 200"p| ), i.e. 
quite comparable to the distance to Mkn 501. 

The second condition for realization of this model - the 
extension of the primary spectrum of Mkn 501 out to 100 
TeV - also is a very robust requirement. Although the 
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Fig. 11. Sketch of the two-stage 7-ray production scenario in Mkn 501. At the first stage, the observed highest energy 
7-rays above 10 TeV with a sharp spectral form ("pile-up") are produced due to comptonization of the ambient optical 
photons of energy ~ 1 — 2 eV by the cold ultrarelativistic jet with bulk motion Lorentz factor F ~ 10 7 — 10 8 . Low 
energy 7-rays, E y < 10 TeV can be also produced due to comptonization of far infrared photons with broad spectral 
distribution. At the second stage, relatively low energy (E 7 < 10 TeV) 7-rays are produced in the blobs moving with 
Lorentz factor of about r ~ 10 in accordance with the standard SSC or external IC models. 



discussion of possible mechanisms of production of such 
energetic photons is beyond the framework of this paper, 
we note that the condition of E max 3> 100 TeV cannot 
be easily accommodated by the conventional (especially, 
leptonic ) models of high energy radiation of Mkn 501. 

Apparently, the "intergalactic cascade" hypothesis 
which simply postulates very low intergalactic magnetic 
fields and very energetic primary 7-radiation of the source, 
cannot provide conclusive predictions for the time corre- 
lations of TeV 7-rays with nonthermal radiation compo- 
nents at other energy bands. At the same time, this hy- 
pothesis allows robust calculations of spectral characteris- 
tics of 7-radiation at different depths of the cascade devel- 
opment. Therefore the most straightforward (and perhaps 
relatively easy) inspection of this hypothesis would be the 
study of spectral characteristics of TeV 7-rays detected 
from sources with different redshifts z. 



5.2. Inverse Compton radiation of ultrarelativistic cold 
outflow 

Although the assumption about the existence of a cold ul- 
trarelativistic outflow with an extreme Lorentz factor of 
bulk motion T ~ 3.5 x 10 7 is somewhat unusual and per- 
haps even provocative (at least it has not been discussed 
in the literature before), it cannot be a priori ruled out. 
A similar scenario most probably takes place, although 
on significantly smaller scales, in environments of pulsars. 
The rotation-powered pulsars eject plasma in the form of 
relativistic winds which carry off bulk of the rotational en- 
ergy. At a distance d < 1 pc the wind is terminated by a 
strong standing reverse shock which accelerates particles 
and randomizes their pitch angles. This results in forma- 
tion of strong synchrotron and inverse Compton nebulae. 
On the other hand, it is generally believed that the region 
between the pulsar magnetosphere and the shock is invis- 
ible because the electrons move together with magnetic 
field and thus do not emit synchrotron radiation. However, 
recently it has been argued that such winds could be di- 



rectly observed through their inverse Compton emission, 
the low-frequency seed photons for comptonization being 
provided by the neutron star in the case of radiopulsars 



like Crab or Vela (Bogovalov & Aharonian 200C) or by 
the optical companion star in the case of binary pulsars 
( Ball fc Kirk 200C| ). If such a cold ultrarelativistic coni- 
cal wind can indeed be produced in the proximity of the 
central rotating black hole, because of possible existence 
of dense photons fields in the inner sub-parsec region, it 
could be a very powerful emitter of inverse Compton 7- 
radiation. The later however would not be accompanied 
by noticeable synchrotron radiation. 



The possibility to disentangle the multi-TeV emission 
with a characteristic sharp pile-up at the very end of the 
spectrum, E 7 ~ m c c 2 L, associated with the unshocked 

the 



jet, from the sub-10 TeV emission associated with 
shocked structures (e.g. blobs) in the jet (see Fig. 
not only may solve the "IR background - TeV gamma-ray 
crisis" but also would allow more relaxed parameter space 
for interpretation of X-rays and the remaining low energy 
(< 10 TeV) 7-rays within the conventional SSC scenario. 
Consequently, this offers more options for interpretation 
of X-ray/TeV 7-ray correlations both on small (t < sev- 
eral hours) and large (weeks or more) timescales. If the 
overall TeV radiation of Mkn 501 indeed consists of two, 
unshocked and shocked jet radiation components (curves I 
and II in Figs. ||, |l0|), we may expect essentially different 
time behaviors of these radiation components. In particu- 
lar, the "unshocked jet" > 10 TeV radiation should arrive 
earlier than the SSC components of radiation consisting 
of synchrotron X-rays and sub-10 TeV 7-rays. Generally, 
all principal parameters of blobs like the magnetic field B, 
the maximum energy of electrons B max , the radius of the 
blob, etc., which define the spectra and absolute fluxes of 
synchrotron X- rays and IC TeV 7-rays, ma y dramatically 
evolve in time ( Coppi fc Aharonian 1999a ). Therefore we 
should expect strong X/TeV correlations, especially dur- 
ing strong flares of the source. These correlations may be 
realized in quite different forms, i.e. it could vary from flare 
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to flare depending on the time-evolution of specific param- 
eters of blobs. At the same time the IC radiation of the 
unshocked jet significantly depends only on the Lorentz 
factor of the jet T. If the latter during short flares re- 
mains unchanged, we should expect rather stable spectral 
shape of the > 10 TeV emission, despite strong varia- 
tions of the SSC radiation (X- and 7-ray) components. 
However, it should be noticed that because of possible ex- 
ternal broad-band infrared radiation fields, the unshocked 
jet radiation component may significantly contribute to 
the low energy TeV radiation as well (see Sect. ||). If so, 
it would enhance correlations between high and low en- 
ergy bands of TeV radiation, On the other hand, it would 
reduce and make more complicated correlations between 
the low energy TeV 7-rays and synchrotron X-rays. 

On larger time-scales, e.g. in a low state of the source 
which could last weeks or months, both the particle den- 
sity and the Lorentz factor of the unshocked jet may be 
smaller than in high states. This would not only reduce 
the luminosity of all components of nonthermal radia- 
tion, but also would make softer the overall 7-ray spec- 
trum. Apparently, the pile-up associated with the radi- 
ation of cold wind would be shifted towards lower ener- 
gies. Although its intensity could be significantly weak- 
ened, nevertheless the search for a relevant spectral fea- 
ture seems to be an interesting opportunity, especially if 
we take into account that at these energies the intergalac- 
tic absorption effect also becomes weaker. Because of low 
photon statistics such studies can be effectively performed 
only by next generation ground-based detectors. 

The straightforward proof of the suggested model 
would be direct detection of TeV 7-ray spectra with char- 
acteristic sharp high energy pile-up. Apparently this re- 
quires nearby TeV sources with d <C 100 Mpc, thus the 
main fraction of > 10 TeV 7-rays would arrive with- 
out significant intergalactic absorption. Unfortunately, all 
currently known BL Lacs are located beyond 100 Mpc. 
This significantly limits the chances for direct detection 
of sharp pile-ups in the spectra of BL Lacs or other type 
of blazars, unless the initial Lorentz factors of ultrarela- 
tivistic outflows in some objects do not exceed L < 10 7 . 

Important tests of the suggested two stage (pre-shock 
plus post shock) scenario of the TeV radiation of jets can be 
provided also by the search for correlations (or lack of such 
correlations) of high energy (multi-TeV) 7-radiation with 
both the low energy (e.g. 1-3 TeV) 7-rays and synchrotron 
X-rays. The low statistics of (heavily absorbed) 7-rays 
above 10 TeV makes the search for such correlations rather 
difficult, and requires ground-based instruments with very 
large, 3> 10 5 m 2 detection areas in this energy domain. 
The new generation imaging Cherenkov telescope arrays 
like CANGAROO-3, H.E.S.S. and VERITAS should be 
able, hopefully, to perform such correlation studies. 
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